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Abstract 



<^ We study asymptotics of the partition function Z^ of a Laguerre-type random 

^^ matrix model when the matrix order A'^ tends to infinity. By using the Deift-Zhou 

r^ steepest descent method for Riemann-Hilbert problems, we obtain an asymptotic 

Cd expansion of log Zj\[ in powers of A^^^. 
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1 Introduction and statement of results 

In this paper, we are interested in asymptotics of the partition function of a Laguerre- 
type random matrix model as follows 

■■■ n (A,-Afc)%"e-^^f-^'(^^)rfAi---rfAAr, «>-l, (1.1) 

l<j<fc<Ar 

where t := (ti,--- ,ti,) and Vt{X) is a polynomial of degree u with positive leading 
coefficient 

V 

Vt{X) := V{X; h,--- ,U) = X + ^ tfcA^ U > 0. (1.2) 

k=l 

It is well-known that partition functions are closely related to orthogonal polynomials. 
Let Pn{x', N, t) be the n-th order orthonormal polynomial with respect to the weight 

Wn{x) := a;"e-^^*(") x E (0,oo), (1.3) 

that is, 

Pm{x; N, t)p„(x; iV, t)a;°e-^^'(")da; = 6m,n, Pn{x; N, t) = 7^ '^^x" + ■ • ■ , (1.4) 



with the leading coefficient 7n > 0. (For the sake of brevity, we shall suppress the N 
and t dependence when there is no confusion.) Then, ZN(t) can be rewritten as 



N-l 



ZN{t) = N\ll{jrY'; (1.5) 



k=l 



see [11]. With the definition in (1.1), it is easily seen that the logarithmic derivative of 



Ziv(t) with respect to the parameter ti is given by 



Q 1 fOO POO ^' 

^log(Z^)=^/ ■■■/ (-iV^Aj.) n (A,-A.)%"e-^^.-^^(^^)dAi---dA;v 

I N Jo Jo j^^ l<j<k<N 

^ /I \ 

=E{-N Y^ Aj) = -A^^E ( ^TrM^ J , (1.6) 

where E denotes the expectation value with respect to the probability measure dfit 

diit = J- (det Mf exp {- NTT[Vt{M)]) dM. (1.7) 

Zn 

In the above formula, dM is Lebesgue measure on N x N positive definite Hermitian 
matrices and Zjsi is the normalization constant such that d/it is a probability measure. 



Like (1.5), the logarithmic derivative can also be put into a form related to orthogonal 



polynomials as follows 

— \ogZr, = -N'J X'p^i\X)dX, (1.8) 

where Pj^ (X) is the so-called "one-point correlation function" 

N-l 

p«(A) = p«(A;t) := lA"e-^^^W5^p.(A)^. (1.9) 

fc=0 



By the fundamental theorem of calculus, we have from (1.8) 

ft ^oo 



ZAr(t) = ZAr(0)exp 



-N^ / / pSi^(A)VVr(A)rfA-rfr 
Jo Jo 



;i.io) 



where Zn{0) is related to the classical Laguerre polynomials and given explicitly below 



N 



Z^iO) = iV-^(^+") H r(j + l)r(j + a); (1.11) 

i=i 

see [H p.321]. 

In the literature, a lot of researchers are interested in asymptotics of partition func- 
tions due to their importance in mathematical physics. For example, for the one-cut reg- 
ular case in the one-matrix model, Ercolani and McLaughlin in [8] applied the Riemann- 
Hilbert techniques to prove that the logarithmic of the partition function has an asymp- 
totic expansion in powers of l/N"^. Later, Bleher and Its |1] used another method to 
obtain this result. Their proof is mainly based on asymptotics of recurrence coefficients 
for the corresponding orthogonal polynomials. For multi-cut case in the one-matrix 
model, only formal asymptotic expansions of the partition functions are derived; see 
[21 [To]. The rigorous mathematical proof is still unknown. For the one-matrix model, 
some people are also interested in cases when there exist some singularities in the model. 
For example, Krasovsky studied power-like (Fisher-Hartwig) singularities in [13]. More- 
over, the asymptotics for partition functions in multi-matrix models are studied, too. 
For instance, one may refer to a series of papers [H [HI [12] done by Guionnet and her 
colleagues. 

In this work, we plan to derive the asymptotic expansion of logZi\f{t) as A^ — )■ oo. 



Based on (1.10), it is sufficient to derive the asymptotics of the one-point correlation 
function pjv (A) for A G [0, oo). As pjv (A) involves orthogonal polynomials Pfc(A), it is 
natural to study the asymptotics of Pk{X) first. This can be achieved by using the well- 
known Deift-Zhou steepest descent method for Riemann-Hilbert problems introduced by 



Deift et. al. in [6l[7]; see also [5]. Note that the idea of using Riemann-Hilbert techniques 
was first adopted by Ercolani and McLaughhn in |8] when they studied an Hermite-type 



random matrix model. Although there are quite a few parameters tk in (1.1), we restrict 



ourselves to the case when the parameter vector t belongs to the following set T(T, 7) 



u-l 



T(r, 7) = U gM'^ : |t| <r,t^ >7^|tj| I for any given r > 0, 7 > 0. (1.12) 

In fact, we will choose T small enough and 7 large enough such that only the one-cut 
case needs to be studied. 

The following is our main result. 



Theorem 1. Assume a = in (1.3). There exist T > and 7 > 0, such that for 



t G T(T, 7), we have the following asymptotic expansion 



fc=0 



where Zn{0) is given in (1.11) and ej{t) is an analytic function of the vector t in a 
neighborhood of for every j. 

Remark 1. It is natural to ask why we only consider such a simple case a = instead of 
general a. The reason is that, when a 7^ 0, nice symmetric properties of the asymptotic 
expansion of pV(A) are lost. Of course, we can still derive the asymptotic expansion of 
log ( Y(^ ) ; but this expansion is given in powers of 1/A^ instead of l/N"^ in the above 
theorem. For more detailed explanation, one may refer to Remark [2] following Lemma [3j 
As mentioned earlier, to prove our main result, we will derive the uniform asymptotic 
expansion of pj^ (A) for A G [0, 00) first. Actually, with its uniform asymptotic expansion, 
we can get the following more general result. 

Theorem 2. Let 0(A) be a C°^ -smooth function and grow no faster than a polynomial 
for A — )■ 00. When a = 0, there exist T > and 7 > 0, such that for all t G T(T, 7), 
the following expansion holds true 

/ e(A)pii)(A;£)ciA~5^-^e,, as N-^ 00, (1.14) 

■^0 fc=0 

where the coefficients Bj depend analytically on t for t G T(T, 7). 

This paper is organized as follows. In Section |2} we present the Riemann-Hilbert 



problem for the orthogonal polynomial Pn{z) in (1.4) and calculate the equilibrium 
measure. In Section |3| we give a sketch about the Deift- Zhou steepest descent analysis. 
Based on the uniform asymptotic expansion obtained for pn{z), we obtain the asymptotic 
expansion for p}^ (z) and show some of its nice properties in Section 4 Finally, in Section 
[5| we prove Theorem [2| which yields Theorem [T] as a direct result. 



2 Riemann-Hilbert problems and the equilibrium mea- 
sure 



.(1)^ 



To obtain the asymptotics for p)^ (X) in (1.9), it is helpful to put it into the following 



a (see [16 



form by using the Christoffel-Darboux formu 

p!i^(A) = ^A"e-^^'W [p'nWpn-i{X) - Pn{X)p'j,^^{X)] ^. (2.1) 



{N 
Tn 



Then, it is obvious that the asymptotics of p)^ {X) is determined by the asymptotics of 
Pn{X) as N —^ oo. To derive the asymptotic expansion oi pn{X), we apply the Deift-Zhou 
steepest descent method for Riemann-Hilbert problems. 

2.1 Riemann-Hilbert problems 

Consider a 2 x 2 Riemann-Hilbert (RH) problem as follows: 

(Y,,) r : C \ M ^ C2^2 jg analytic for C \ [0, cx)); 

(Yfo) Y{z) possesses continuous boundary values on (0, oo). Let y+(x) and Y^{x) denote 
the limiting value of Y{z) as z tends to x from above and below, respectively. They 
satisfy 

/ 1 ^a -NVtix) \ 

y+(x) = F_(x) ^ ^ forxG(0,oo); 



^ 

(YJ for ^gC\[0,cx)), 

Y(z) 



~z 



z"" 
z 



'N 



as ^ — )■ oo; (2.2) 



(v.,K,,^o (;;;;;(;)) 



as z ^ 0,z E C\[0, oo), where ria{z) is defined by 



1, if a > 0, 

Va{z)={\og{l/\z\), ifa = 0, (2.3) 

|z|", if - 1 < a < 0. 

According to the significant results of Fokas, Its and Kitaev [H] , the solution of the 
above RH problem is given in terms of the monic polynomials 7r7v(a;) orthogonal with 
respect to w{x). This establishes an important relation between orthogonal polynomials 
and Riemann-Hilbert problems. 



Lemma 1. (Fokas, Its and Kitaev 19]). The unique solution to the above RH problem 
for Y is given by 

Yfz) = ( ''''^^^ C[7Tnw]{z) \ 



where w is the weight function given in (1.3), cat = —2Tii'^]^_^ and 



C[f]iz) := ^ r P^dC zeC\[0,oo), 
^TT? Jo Q- z 

is the Cauchy transform of f . 

2.2 Equilibrium measure 

The equilibrium measure plays an important role in the Riemann-Hilbert analysis and 
we wish to calculate it explicitly. Recall that the equilibrium measure fiy is the unique 
minimizer of the following weighted energy 



-^v(/^) = / / log": 7dfx{x)dfx{y) + / V{x)dfx{x) (2.5) 

J J \x-y\ J 



among all probability measures on [0, oo). For our problem, we have the explicit formula 
as follows. 

Theorem 3. There are % > and 7o > such that for all < T < Tq and 7 > 70, the 
following holds true. If t E T(T, 7), then we have 



1 (3 — X 
d[iv = i^v{x)dx, ^pv{x) = —X{o,i3){x)d h{x), (2.6) 

where h{x) is a polynomials of degree v — \. This polynomial is strictly positive on the 
interval [0, (5] and defined as 



^w-s^i.S"''*^'?^- ^°'- ^^c^i"-"'- p'^) 



where Vt{y) is given in (1.2), V^ is a positively oriented contour containing [0,/3] and z 



in its interior. The endpoint (5 is determined by the equation 

V;ix)J-^dx = 27r. (2.8) 

y p — X 

Readers may compare the above theorem with Theorem 3.1 in [8] and find some 



similarities. The main difference is that equilibrium measure in (2.6) is supported on 



[0, /3] and possesses a square root singularity at 0, while the measure in [8| Thm 3.1] is 



supported on [a, /3] with a < < (3 and vanishes hke a square root at both endpoints. In 
fact, our proof of the above theorem is based on [8l Thm 3.1] and a nice relation between 
the equihbrium measure /iy and the other one uy , which is the unique minimizer of the 
weighted energy 

Iy^^,y,(u) = J J\ogj^duix)duiy) + J^^duix) (2.9) 

among all probability measures on M. The relation is given in the following lemma. 
Lemma 2. Let ipy be the density of fiy CLnd (pv be the density of vy. We have 

(j)y[x) = \x\tljy{x^) forx e R. (2.10) 

Proof. See [3l Lemma 2.2]. D 

With the above Lemma, we can prove Theorem |3] 
Proof of Theorem^ Let Wt{x) be the following potential supported on the whole real 



axis 



Wt{x) = \vt{x^) = \x^ + Y,%^^'' ^eM. (2.11) 



2 ' ' 2 ^^2 

fc=i 



Then according to [H Thm. 3.1], its corresponding equilibrium measure is given by 

dvw = 7^Xia,^]{x)\/ix - a){(3 - x) h{x)dx, (2.12) 



27r 
where h{z) is a polynomial of degree 2z/ — 2 and defined as 



"«^ii^-^m;i^- 



^{s-<i){s-P) 



Here the integral is taken on a positively oriented contour F^ containing (a, /3) and z in 
its interior, the endpoints a and /3 are determined by the following two integrals 

"-'<-' ^- n ,„d t^m^is = 2.. (2.14) 



(s-a)W-s) ■>' J{s-a)W-s) 



Recall (2.11), Wl{x) is an odd function. Then the above formulas gives us a = —/3 with 
/3 > and 

di^w = ^X[-/3,^](a;) V /3^ - x^ h{x)dx. (2.15) 



Using (2.10), we obtain (2.6) with/3 = /3^ and h{x) = h{^/x). In fact, one can easily verify 
that (2.8) and (2.14) are consistent with the relations (2.11). To see that h{x) = h{^/x) 



is a polynomial in x of degree i^ — 1, we recall that W'^^{s) is a polynomial of odd powers. 



Then calculating residue at oo in (2.13), it is easily seen that h{z) is a polynomial of 



even powers. As a consequence, h{x) 
finishes the proof of our theorem. 



h{^/x) is a polynomial in x of degree z/ — 1. This 

D 



3 Deift-Zhou steepest descent analysis 

In the standard Deift-Zhou steepest descent analysis, one introduces a sequence of trans- 
formations: 

Y ^U ^T, 

such that T{z) satisfies a RH problem with simplified jump conditions. Then, when N 
is large, some parametricies T'^'^^z) are constructed to approximate T{z) in different 
regions of the complex- 2; plane. As the above transformations are revertible, one gets 
large- A^ asymptotics for Y{z) from T^'^^z). 



For our problem, since the equilibrium measure in (2.6) is supported on one interval 



this is the so-called one-cut case. The Riemann-Hilbert analysis is similar to that done 
by Vanlessen in [TTj. So we only give a sketch for the completeness of this paper. The 
interested readers may refer to pT| for details. 

Normalization: F — t- f/. The first transformation in the Deift-Zhou steepest descent 



analysis is to normalize the large-z behavior of Y[z) in (2.2) and make it tend to the 



identity matrix. To achieve it, we introduce the following ^f- function 
log(2: — x)il)v{x)dx 



9[z) 



1 
2^ 



/3 — X 

\og(z — x)\l h(x)dx, z e C\(— 00, /5], 

V X 

(3.1) 

where the principal branch of the logarithm is taken. It is obvious that g{z) behaves like 

log z when z is large. Then, the first transformation Y -^ U is defined as 



U{z) 



'-''^Y(z)e~^'^^^^''^e- 



zeC\]\ 



(3.2) 



such that the large- 2; behavior of U{z) is U{z) = I + 0{z ^) as z -^ 00. Here a^ is the 
'1 



Pauli matrix 



-1 



Opening of the lens: U ^f T. In the second transformation U ^ T, we deform the 
original interval [0, 00) and open lens. With this transformation, the rapidly oscillatory 



jump matrices for U will be reduced to jump matrices who tend to I exponentially except 
in the neighbourhood of (0,/3). 

Before we introduce the transformation, we need some auxiliary functions. Define 



^W := ^V ^^W' for z e C\[0,/3], 
where the principal branch of the square root is taken, and 



(3.3) 



^{z) := - / ^{y)dy = -- 



y 



h{y)dy for 2; G C\(— oo, /3], 



(3.4) 



where the path of integration does not cross the real axis and the principal branch of 
the square root is taken. The ^-function satisfies the following properties. 



Proposition 1. For x G M, we have 

^+(x) — ^_(x) = 27ri /or X G (—00,0), 

2e+(x) = -2e-(x) = g^{x)-g4x) for x e (0,/3), 

^{x) < for X e {l3, 00). 

Moreover, there exists a 6 > such that 

Re^{z)>0 /or < llm^l < 5, <Rez < jS. 

and 

2az) = 2g{z) - Vt{z) - ly, forzeC\ [/3, 00), 

where ly is a constant. 

Proof. The proof is similar to the analysis in ^71 Sec. 3.4]. 
The second transformation [/ — )■ T is defined as 

[/(z), for z outside the lens-shaped region 

U{z) \ ^^^^^ ^ I , for z in the upper lens region. 



(3.5) 
(3.6) 
(3.7) 



(3.8) 
(3.9) 

D 



nz) := 



U{z) 



1 

1 0' 



(3.10) 



for z in the lower lens region; 



see Fig. [1} Then, T satisfies a RH problem with the following jump conditions Jt{ 




Figure 1: Contour S^ = |Ji=i ^i 

^ G S1US3 

/ z°'\ 

Jt(z) = ({ 1, ^e(0,/3) 

2;e (/3,oo). 
Of course T{z) is still normalized at 00, that is T{z) = 1 + (z^^) as 2; — )■ 00. 



1 

z° 

-Z-" 

,0 1 



(3.11) 



3.1 Outside parametrix 

According to the properties of ^(z) in Proposition [11 we can choose Si and S3 such 
that Re^(z) > for 2 G Si U S3, and Re^(a;) < for x G (/?, 00). Thus, Jt{z) -)■ / 
exponentially as A^ — )■ 00 for ^ bounded away from [0, /?]. It suggests that, for large A^, 
the solution of the RH problem for T{z) behaves asymptotically like the solution of the 
following RH problem for T^'^^z): 

(a) T(°°)(z) is analytic in C \ [0,/3]; 



(b) forxG (0,/3), 



T^r\x) 



t[^Hx) 




—x~ 







(3.12) 



(c) T(°°)(z) = / + 0{z-^), as z^oo. 
The solution to the above RH problem is given explicitly as follows 



2^(00), 



a{z)+a ■'-(2) 

^ ' a{z)-a-^z) 

-2i 



a{z)—a ^(z) 
a{z)+a-'-(z) I -^C^^ 



ZGC\[0,/3] 



(3.13) 



10 



with 



and 



a{z) 



P) 



ioT zeC\[0,(3] 



z^ 



D(z) 



Z2 



[z). 



[z) = 2z - (3 + 2z^z - (3)^ ioT z e C\[0, 13], 



(3.14) 



(3.15) 



where principal branches of powers are taken in the above formulas. 

Note that the above parametrix is not valid in the neighbourhoods of endpoints and 
P, since the jump matrices do not tend to / as A^ — )■ oo. So, some local parametrix con- 
structions are needed near there two endpoints. As the equilibrium measure has different 
behaviors near and (3 (see Theorem [3]), different parametrices appear in their neigh- 
borhoods. More precisely, we will have an Airy-type parametrix near /3 since the density 
of equilibrium measure vanishes like a square root at (3; and a Bessel-type parametrix 
near as the density has a square root singularity there. 



3.2 Local parametrix near (3 

This parametrix is constructed in terms of Airy functions as follows 

T^^\z) = E('^)(z)^(^)(/(z))e-^«(^)"3^-^' 
where 



jao-a 



E('^)f2) = T(°")f2)25°-3ef-3 



( 



for z G C\(— oo,/3], and 



3A^ 



(p{s)ds 



^(/3)( 



27re 12 X < 




1 -1 



V2V 1 



/(^; 



4 




h(s)ds 




e e 



e 6 



0-3 



0-3 



e 6 




-0-3 



zel 
zell 
zelll 
zeW 



(3.16) 

(3.17) 

(3.18) 



(3.19) 



with CO = es'^*; see [T71 eq.(3.62)]. Here the regions I— IV are depicted in Fig. \2\ For 
properties of Airy functions, one may refer to |15j . 



11 




Figure 2: Regions for the Airy parametrix with the angle a G (f , vr) 



3.3 Local parametrix near 

This parametrix is constructed in terms of Bessel functions as follows 

where 



E(o)(z) = (-l)^T(°°)(2)(-2)l"-=^^ /"^ 






for 2; G C\[0, 00), and 



iV 



ip{s)ds 



N 



S-/3 



h{s)ds 



^(o)(, 



-27ri25/;(225) -2^5ir;(2^i) 
-7rzi(/7i^))'(2(-^i)) 7rzi(ffFy(2(-zi)) 



zeV 

e—"'-", z G ir 
e-— ""^^ z G III'; 



(3.20) 

(3.21) 
(3.22) 



(3.23) 



see [T71 eq.(3.81)]. The regions I'— III' are depicted in Fig. 13} Here /„ and Ka are modified 



r(l) 



r{2) 



Bessel functions, and Ha and Ha are Hankel functions of the first and second kind, 
respectively. For properties of these functions, one may refer to [15] again. 



12 




Figure 3: Regions for the Bessel parametrix with the angle a G (0, vr) 

3.4 The error 

Now we have the following approximation for T{z) 

't^^\z) for^Gf/A^T, 
T(^)(2) = <J TW(z) for z G f/ASr, 
T(°^)(z) elsewhere, 



(3.24) 



where U^ and U^ are two disks with constant radius 5, centered at (3 and 0, respectively. 
To see the difference between T{z) and T^'^^z)^ we define 



S{z) ■.= T{z)T^^\z)-\ 
Then, S{z) satisfies a RH problem as follows. 

{So) 



(3.25) 



S{z) is analytic in C\Ss; see Fig. |4 
S^{z) = S^{z)Js{z), where 



T(o)(2)T(°°)(^)-i 
Js(z) = ( T^P){z)T^-^\z)-^ 

T(^\z)Jt{z)T<^'^\z)-' for z G Ei,5 U ^3,5 U ^4,^; 



for z G dUs 
for z G dUs 



(3.26) 



(5c) S{z) = l + 0{1) asz^oo. 

Based on the formulas for Jt{z) and T^°°\z) in (3.11) and (3.13), one can see Js{z) 



is exponentially close to the identity matrix for z G T^sXidUs U dUs) as A^ — )■ oo. For 
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Figure 4: Contour E^ 



z G dU^ U dUi, Js{z) possess the following asymptotic expansions as A^ — )■ oo 

T(«(z)T(-)(^)-i ~ / + 5^ ^-/ri^), for z G dUs, 



J,(z) = < 



k=l 

oo 



T(0)(^)T(-)(z)-i ~ / + 5^ ^4°^ W, for z G 9^,, 



(3.27) 



fc=i 



where Jp, (z) and J^ (z) are meromorphic functions in Us and f/^ with a pole at /3 and 
0, respectively. Moreover, they have explicit expressions as follows 



jr(^) 



1 



2ij;^is)dsr 



rp{oo). 



Z)Z^ 



ca 



{-lY{Sk+tk) {Sk-tk)i 

{-lY+^{sk-tk)i Sk + tk 



^-f<x3y{oo)(^^-l 



(3.28) 



and 



J, 



(0) 



(^) = fe^^44#r'-'(.)(-.)^ 



'(-I)'' ^^,2 I 1 



^0-3 



m;vis)d 



{a' + lk-D (fc-l) 



(-i^M^-i)^ H«^ + ^^-i), 



X (-z) 2 



-32-(oo). x-l 



(3.29) 



where 



Sk 



r(3A;+|; 
54'=A;!r(A; + 



-)' 



6fc + 1 , , , (4a2 - 1) ■ • • [4a2 _ (2A; - l)^] 
tk = --, 7Sfc> {a,k)=^- ^ ^ ^ ^J 



6A;-1 



22fc^! 



for A; > 1; see [T71 eq.(3.75)&(3.97)]. Because the jump matrix tends to I as N -^ 00, 
by the standard analysis in [7], the RH problem for 5* has a unique solution when N is 
sufficiently large. Moreover, S{z) satisfies the following asymptotic expansion 



S{z)r^ I + ^N-^Sk{z), asN^oo 



(3.30) 



k=l 
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where Sk{z) are bounded functions which are analytic in C \ {dUs U dUs) and can be 
computed exphcitly. 

There is an important observation about exphcit forms of Sk{z). From formulas for 
T(°°)(z), Ji^\z) and jf^(z) in ( [3l3| , ( p^ and ( |3^ , one can see that they satisfy 
nice symmetric properties when a = 0. Due to these properties, we have the following 
results for Sk{z), which are crucial in our subsequent analysis. 

Lemma 3. When a = 0, we have 

sl {z)I + s\ {z)a2, when k is even, 
s], {z)a'j, + s], (-2)o"i, when k is odd, 

where Ui and 02 are the Pauli matrices 

'0 i\ /o -r 

Proof. When a = 0, the outside parametrix T'^^^z) in (3.13) is simply given by 



Skiz) 



0-1 



(3.31) 



(3.32) 



rrioo), . ^ a{z) + a \z) a{z)-a \z) 
^ ' 2 2 



0-2 • 



(3.33) 



m. 



Then, for J^ (z) in (3.28), we have 
2fc-i 



r(/3) 



jn^) 



u:v{s)dsf 



rp{oo). 



{sk + tk)I - {sk -tk)a2 



T^'^Hz) 



■)k-l 



Ub ^is)ds 



[Sk +tk)I - {Sk -tk)(T2 



(3.34) 



and 



Ji%) 



)fc-i 



u;v{s)dsY 
{!;^{s)dsY 



rpico)^. 



{Sk + tk)cr3 + i{sk-tk)cri 



when k is even, 



^(oo)(^)-l 



~(a\z)tk + 



Sk 



a^{z) 



0-3 - a {z)tk - 



Sk 



a^{z) 



lai 



, when k is odd. 
(3.35) 



Similarly, for J]. '{z) in (3.29), we get 



J^^\z^ (0,fc-l)(-l) 



k r 



and 



J^\z) 



mo^is)ds) 



(o,fc-i)(-i) 

2^(/X.)rf.)' 



1 1 

2 ~ 4A? 



k-^\^2 



when k is even, (3.36) 



{2k-lfa\z) {2k-l){2k + l] 



8k 



8ka^(z) 



0-3 



■ (2A:-l)V(z) (2A:-l)(2fc + l) 

+ I ^^1 \ —r 7T^-^ I 'iCTl 



8k 



]ka^{z) 



(3.37) 



when k is odd. 



15 



According to the above formulas, the jump matrix Js{z) in (3.27) always satisfies an 
asymptotic expansion with symmetric properties. That is, its even terms are given as a 
linear combination of / and o"2, while its odd terms are expressed as a linear combination 
of (Ti and (73. Therefore, using the same idea as in the proof of [HI Lemma 3.5], it is not 
difficult to prove our results by mathematical induction. This completes the proof of our 
lemma. D 

Remark 2. When a 7^ 0, it is still possible to calculate J^ \z) and J^ (z). But now the 
formulas are more complicated. For example, when k is even, we have 



Ji'\^) 



{j;^{s)dsf 



u;^{s)dsf 



,a + a ^ ^2 f 



,a — a 



-1 



2 / -iz'^'D^A " \ a- -a 



a^ - a-2 fz"D-^ 




z-'^D^ 



an2/i-a\ ^2 

+ 





-2 ^z-"D^ 

2"L)-2 







where a and D denote the functions a{z) and D(z) in (3.14) and (3.15) for brevity. 
For such kind of JI (z) and j|, (z), the pattern for Sk{z) in (3.31) is no longer valid. 



Therefore, it is likely that one can only obtain the asymptotic expansion in powers of 
1/A^ instead of l/N'^ in Theorem pi This is also the reason why we focus on the case 
a = in this paper. 

4 Asymptotic expansion of p\J{z) 

With the Deift-Zhou steepest descent analysis done in the previous section, we are ready 
to obtain the asymptotic expansion for p)^ \z). Reverting the transformations in (3.2), 



(3.10) and (3.25), we have 

Yiz 



e 2 



-0-3 



U{z)e 



N(giz)-'-y-)a3 



e^^-'Tiz) 



r{z) 



0^ -,N(g{z)-'f)a3 



'^"'S{z)T^^\z) 



r{z) 



^N{9(z)-'-f)a^ 



(4.1) 



where 



r{z) 



0, for z outside the lens-shaped region in Fig. [T] 

-|-g-2Af^(2;)^ for z in the upper/lower lens region in Fig. [I| 
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Recalling (2.1) and (2.4), we get 

(1), . e-^^*(^) 



p'N'iz) 



2niN 



[Y|,iz)Y,,iz)-Y^,iz)Y^,iz)]. 



(4.2) 



For z in the neighborhood of /3, since Y{z) involves the same Airy-type parametrix, we 
have the same formula as that in [8; eq.(4.4)]. For z in the neighborhood of 0, as the 
Bessel-type paramatrix is constructed, the computations are different. First, we get the 
following expression. 

Lemma 4. Let t belong to the set T(T, 7) in Theorem^ With a{z) and f{z) defined in 



(3.14) and (3.22), we have 



pSH^) = I (^ + {j) f"'h{2f'")m^") - ^ 



/o^(2/V2) _ (j;(2/V2)) 



ANni 



(5'(z)i?(°)(z)^(/))n(5(^)i?(°)(^)^(/))2i-(5(^)5W(^)^(/))n(5'(^)5(°)(^)^(/))2i 

(4.3) 



(3.26), 



for z in a fixed-size complex neighborhood of 0. Here S{z) satisfies the RH problem in 

for zeC\ {-00,0], (4.4) 



^ z := 



/o(2^V2) 
-2mz^/H'o{2z'/^) 



-^Ka{2z^l^) 
-2z^'^K'^{2z^l^\ 



and 



EW(^) = f ^^;^^ "'\}'l]f{z)t(2^y 



Proof. According to the parametrix in (3.20), we have from (3.9) and (4.1) 



Y{z) = e^"3^(^)E(°)(^)^(/>))e^^'(")"3^ 



(4.5) 



(4.6) 



where 



E^'\z) ={-lfT 



(oo). 



V2 \i 



f{z)^{2n)^ 



\N 



la 



\N 



, \f{z)^{2n)^ 



-B^'\z). 



V2 \ia a ^ J "■ ' "■ ' 72 

As Yl-^{z) and Y2i{z) appear in the formula for p]^'{z), we need to calculate the derivative 
for (4.6). From the formula for E^'^^z) in the above formula, we have 



d 

dz 



E^^\z) 



-1 



iTV 



V2 



-ia ^a' \ -^ 



la —a (X 



- I ^^ , a la 



1\ / 1 £/ 



la a 



f'f- 



iff-^ 



a' f 
a 4/ 



E^i 



z as. 



(2vr)^ 
(4.7) 
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To calculate -^^(z), we recall that both Ia{z) and Ka{z) satisfy the modified Bessel 
equation 

z^u"{z) + zu'{z) - {z^ + a'^)u{z) = 0. 

Then we get 



Combining the above formulas (4.6)-(4.8), we obtain (4.3). D 



(1)^ 



To derive the full asymptotic expansion for p}^'{z), some more detailed computa- 



tions about the last term in (4.3) are needed. Based on the pattern for the asymptotic 



expansion of S{z) in Lemma pi we get the following result. 
Lemma 5. We have 



ANtTiI'"" ----- ---- -- -v.....v~- -W..21 

~ Y: N-^~a,{z)a^~fyHl{2f'l')+ Y. ^-^•^.(^) ^~'^'^°^y~'^'^' (4.9) 

j even, j>2 j even, j>2 

j odd, j>3 

where aj{z), bj{z) and Cj{z) are analytic functions in both z and t, for z in a fixed-size 
neighborhood of and t belong to the set T(T, 7) in Theorem^ 

Proof. Recalling Lemma |3| we have 

k odd, k>l k even, fe>2 

which gives us 

S\^)- E (4'^'^3 + 4'^VOiV-^ + Y. {sf I + sf a,)N-\ as AT -> 00. 

k odd, fc>l k even, fc>2 

Consider the first column of 

B^^\z)nf) = h' *) , (4.10) 

\m2i *J 
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where B^'^^z) is given in (4.5). It is easily verified that 



1/2N2 



ml, - ml, = 4nf'/' (a'P,{2p/') - ^ 



mnm^i = 2m f'/' | a'P,{2f'/') - &p^ 



(4.11) 



From the above formulas, we have 



]\[2k-l "^ Z^ 



"^21 



k 



hk-imn - S2k-im2i 

J\f2k-1 



k 



N'ik 



i hkmii + 52fc^2i 

JY2fc 



(4.12) 



and 



(55(°)^)n(^'S(°)^) 



21 



mil 



E 



S2fc-iy»ll + t2k-im21 

]\[2k-l 



E 



S2fc?^ii -it2km2i 



J^2k 



X 



Z^ /V2fc-1 ~^ Z^ 



]\[2k 



(4.13) 



Substracting the above two formulas, one can see that all odd terms of 1/A^ are given 
in terms of m\, — m\, and miim,2i- On the other hand, all even terms of 1/A^ are given 



in terms of m\, + m^i- Recalling the properties of mn and mi2 in (4.11), one gets 



(4.9). And analyticity of the coefficients aj{z), bj{z) and Cj{z) follows from the analytic 
property of 5* (2;). D 



5 Proof of the main theorem 

Now we are ready to prove Theorem |2] based on the representations of the one-point 
correlation function Pj^ (z) obtained in the previous section. We would like to carry 
on our asymptotic analysis in the neighborhood of and /3 separately. To achieve this 
object, let us introduce a partition {xo, X13} of unity for (0, 00) as follows 

Xoiz) is C°° with < xo(^) < 1 

supp Xo C {0,z* +e), Xo = 1, for ze[0,z*-e) 

and xisiz) = 1 - Xo(^), 
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where z* is a fixed point in (0,/3). Note that Pj^ (z) is exponentially small for large N 

POO 

when z is bounded away from [0, /3]. Then, instead of considering / Q{X)p)^'{X)dX, it is 

Jo 

enough to study / Q{X)p}/{X)d\ for some small A^-independent 5. Recall that 0(A) 

Jo 
is a C°^-smooth function and grows no faster than a polynomial for A — )■ oo. According 

to the partition of unity above, we rewrite p]y {z) as 



p«(^) = Xo(^)p!i'°H^) + Xp{z)p%^'\z), (5.1) 



and obtain 



e{z)p^^\z)dz= xo{z)e{z)p'^i''\z)dz+ x^izMz)p^i'^\z)dz. (5.2) 

Jo Jz*-e 

Because p)^' (z) satisfies the same Airy-type representation as that in [8^, eq.(4.4)], the 
second term on the right-hand side of the above formula has an asymptotic expansion 
in powers of 1/A^^ according to the results in [8]. Then it is sufficient for us to consider 
only the first term 

\oiz)Qiz)pi^'\z)dz. (5.3) 



According to the asymptotic expansion of p]^ in Lemma k^ and p| there are four types 
of integrals in (5.3) as follows 

(i) j^jf^'eiz) [imy^] - (/^(2/V2))2] f'dz, 

(ii) N-^+' Jf^' 9{z)Io{2fy')ro{2f'^')dz ijodd); 
(iii) N~^+' J^'^' 9{z)Ii{2fy')dz {3 even)- 
(iv) N-^^^g'^'e{z){I',Cl]^l'')fdz {3 even)- 



where Qi^z) is a general infinitely differentiable function of z and is compactly supported 
within (0,2;* +e). 

Next, we will show that all integrals of these four types possess asymptotic expansions 
in powers of l/N"^. Before going to the proof of their expansions, we would like to mention 
the following relations among these four types integrals. 

Proposition 2. // integrals of the types (i) have an asymptotic expansion in powers of 
l/N"^ , then integrals of all other three types have asymptotic expansions in powers of 
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Proof. First, for an even integer j, we have the following relations between type (i) and 
(iv) from integration by parts 



AT-^'+i / e{z){I'^{2f^/^)ydz = N-^+^ 



d{z)_f 

r 



{mr")f~r^fdz 



.N-j+i 



_iv-j-i 



f' 

e{z)T\' 1 



/o^(2/V2) _ (j^(2/V2)). 



dz 



T / T 



Il{2f'/') - {I',{2f'/')f f'dz 



Note that, in the integration by parts, the boundary terms don't appear because 9{z) 
is compactly supported within {0,z* + e). Moreover, f'{z) is analytic and nonvanishing 
throughout the region of integration so that the integral is well defined. Lastly, 

t'{z) := N-'fiz) 

defines an analytic nonvanishing function which is independent of A^; see the definition of 



f{z) in ( |3.22[ ). As the quantity {9{z)t/t')' (1/r') is infinitely different iable on [0,z* + e] 
this suggests type (iv) integral can be reduced to type (i). 
Similarly, for the type (iii), we have 



iV-^+i r ^\{z)ll{2r'^)dz 
Jo 

"^*+" e{z 



._N-J+i 



--N-^-'^ 



/' 







r 



Il{2f'l') - {r,{2f^/')f\ f'dz + AT-^'+i 



9{z){I',{2f^'')fdz 
e{z){r,{2f^'')fdz. 



The two terms in the last equation are integrals of the type (i) and (iv), respectively. 
Finally, suppose that j is an odd integer, then, for the type (ii), we have 

2Ar-^+^ r ^%(^)/o(2/V2)/^(2/V2)rf, = 2N-^^' f ""' ^M^Io{2r/y,{2r/')f~'/'f'dz 
Jo Jo f 

I 



.AT-i+i 



0{^)f 



1/2 



/' 



ii{2f''')dz = -N-^ r ^' 

Jo 



e{z)T 



1/2 



ii{2r")d 



z. 



where the last integral is of type (iii). Thus, there four types of integrals are equivalent 
and our results follow. D 

Now we focus on the type (i) integral and derive its asymptotic expansion. We have 
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Lemma 6. The type (i) integral 

'0 L ^ 

has an asymptotic expansion in powers ofl/N"^. 



1 rz+£ 
N 



(5.4) 



We will apply integration by parts and prove the above Lemma by induction. To 
achieve this goal, we need a few preliminary results. Define 



C^z) 



N 



(3-s 



h{s)ds, 



then, we have from the definition of / in (3.22) 

4/ = e±-X' = -C' and 4/ = -SCV- 



(5.5) 



(5.6) 



With the help of the above formula, (5.4) can be rewritten as 

1£*^%(.) [/o^(2/V^) - (/;(2/V2))^] fdz = -A^^ 
where -Fo(C) is an analytic function in ( defined below 



z*+e 



e{z)Fo{C)Cdz, (5.7) 



(5.^ 



Proposition 3. The function Fq satisfies the following asymptotic expansion as ( -^ oo. 

/ oo oo oo \ 

FoiO - e" Y. ^^^~'' + E ^^^~'' ^i^(2C='/2) + J2 ^^''^'" cos(2C='/2) , (5.9) 



, j=0 



i=l 



i=l 



where Ci, Ci and Ci are constants which can he computed explicitly. 

Proof. It is well-known that the Bessel functions Ju{z) and its derivative satisfy the 
following asymptotic expansion as z — )■ oo 



^o(^) 
where 



TIZ 



1/2 



2 

-KZ 



1/2 



7r\ / TT 

COS \ z — - Wi — sin hz — T ] W2 



sm ( 2; - - j W3 + cos (^z - - ) u;4 ) , 



arg2;| < tt — 5, (5.10) 
I argz| < TT — (5, (5.11) 



(-1) T^' ^2~2^(-i) 



fe=0 



Z^ 



y2k^\ ' 



W-i 



E(-i 



jMk_ 



W4 



E(-i 



A;=0 



fe=0 



fe=0 



,fc ^2fc+l . 
' 22fc+l ' 

(5.12) 
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see [15], §10.17]. The coefficients a^ and bk are known explicitly 

_ (4 - 12)(4 - 32) ■ ■ ■ (4 - {2k - 1)2) _ (-l)'^-i((2A; - 3)!!)2(4A;2 - 1) 

with ao = 1, &o = 1 and bi = 3/8. 



Substitute the expansions of the Bessel functions (5.10) and (5.11) into (5.8), we get 
the expansion for Fo{() as follows 



i^o(C) ~ -C'/'(/r + /f sin(2CV2) + f§^ cos(2CV2)), (5.13) 

7r 

where /g , /g and /^ are asymptotic series valid for C — )• oo, 

/o°^ = Wi + ^2 + W3 + W;4, fP =wj-wl-wl+ WJ, fP = 2{wiW2 - W3W4). 

Now using definitions (5.12) of the asymptotic series {wj}^^;^, we see that /q is an 



(0) 



asymptotic expansion in powers of C~ starting with the constant term. Similarly, fg 
is an asymptotic expansion in powers of C~^ starting with (~^, and /^ is an asymptotic 
expansion in powers of the form C"^''"^''^) J = 0, 1, ■ ■ ■ . This finishes the proof of our 
proposition. D 

We need one more lemma as follows. 

Lemma 7. // a function F{() is C°° on [0, 00), and possess the following asymptotic 

expansion 

00 

^(C) ~ E c^C'^'-''^' trig(2C^/2)^ as (^ 00, (5.14) 

j=0 

in which j G N, and trig(-) denotes either sin(-) or cos(-), then the function J_°° F{()d( 
possesses the following asymptotic expansion as s —t- 00. ■ 

FiOdC ~ 5^c:s-(^+^)/2-3./2 trig(2C^/2). (5.15) 

i=0 

Proof. The proof is similar to that of [H Lemma 5.5]. D 



Now let us derive the ffist two terms of the asymptotic expansion of (5.4), together 
with its error term. 

Proposition 4. We have 

~IN f^^ ^(^)^o(C)C'rf^ = eo + ^ei + A,, (5.16) 

where Cq one? ei are independent of N and \A2\ < CN^'^^^ for a constant C . 
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Proof. We rewrite the right-hand side of (5.7) as 

3 '"*+' 



AN 



Jo 



where 



and 



eo = - 



4N 



co9{z)C/'Cdz, 



Ai = -^ £ ^^e{z) (Fo(C) - coC/') Cdz. 



(5.17) 

(5.18) 
(5.19) 



Ahhough there is a factor A^ in the leading term eo, it is actually A^-independent 

3 '"'+' 

Co 



eo 



I' "" c,e{z)C''C'dz = -^ y ^' c,e{z)de/' 



4 





9{z)\ -h{z)dz- 







z 



(5.20) 



see the definition of C,{z) in (5.5). To get the next term, we apply integration by parts 



in (5.19) and obtain 

3 '•"*+" 



m 



£ "" e{z)F,{OCdz = eo + ^iy^ e\z)F, 



{Qdz 



(5.21) 



with 



F[{() = F,{C)-c,e". (5.22) 

We don't have any contributions from the boundary terms because 9{z) is compactly 
supported within (0, z* + e). Next, we repeat integration by parts twice to produce the 
higher order terms 



where 



and 



A, 



3 

AN 



ei 



m 



[ e{z)Fo{OCdz = eo + eiAr-2 + A^, 
Jo 



AN 



1 /I 



c vc 



}■ (?"'<-■ 



:e'(z] 



zfe^'dz, 



/-I ' 



{FsiO - cfe/')dz. 



(5.23) 



(5.24) 



(5.25) 



Here the constant c\ in (5.24) and (5.25) is the leading coefficient of F^^^Q defined later. 
This function -F3(C) satisfies the following relation 



F^(C) = Fi(C), F^(C)=F2(C). 



(5.26) 
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Again, since ( is of order N"^/^ (cf. (5.5)), ei is also independent of A^. To get the 
approximation of A2, we recall the relations among Fi{() in (5.22) and (5.26). Define 
Fi{() as follows 

/■oo 

(5.27) 



/oo 
(Fois) - cos'/')ds. 



Since the integrand in the above integral behaves like 0(C~^''^) as C ^ 00, then -Fi(C) 
exists and satisfies (5.22). Moreover, according to Lemma [TI it satisfies the following 
asymptotic expansion as C — )■ 00 

00 00 00 

FiiO - Y. <^Cr"'-'^-'^ + E <^Cr"'-'''~'^ sin(2C^/2) + ^ ^C'' cos(2C^/2), 

(5.28) 
where c{ = — 2ci/3. Similarly as (5.27), we define 

/oo /"OO 

Fi(s)ds and F^{Q := (F^is) - cf\s-'/^)ds - 2cf\'/^ . (5.29) 

These two functions exist and satisfy the following expansions 



^2(C)~E^S\ 



(2)^-l/2-3(i-l) 



+ E^f^C-^/^-^^sin(2C^/^) + 5^£fC 



~(2)^-2-3(i-l) 



cos(2C 



3/2N 



i=l 



(5.30) 



i^3(c)-E^?^c^'"''"^ + E^f^c- 



5/2-3(J-l) 



sin(2C^/2) + J2 ^k''^'' cos(2C3/2)_ 



i=l 



i=l 



i=l 



(5.31) 

As F-i{() — c\ C^/^ is uniformly bounded, we get the \A2\ < CN^'^^'^, which completes 
the proof. D 

With the above preparations, we are ready for the proof of Lemma |6| 

Proof of Lemma \^ We prove this result by induction. For a positive integer k, define 
F3fc+i(C), -^3/fc+2(C) and -^3^+3(0 as follows: 



(Fskis) - cf^h^'^)ds, F3fc+2(C) = - / F^k+i{s)ds (5.32) 



and 



-^3A:+3(C) 



{F,,^,{s)-cf^'\s-''')ds-2cf'^'\''\ 



(5.33) 



Assume these functions satisfy the following asymptotic expansions for all integers k < 
j — 1 as C — ^ 00: 

00 
i^3.+i(C) - e" E <t^''e~'' + /f ^'^ sin(2C='/2) + f§'^'^ cos(2C^/2), (5.34) 
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' ^(3fe+l) _ >-(3fc+6)/2 y^oo ^(3fe+l) >-3(i-i) 

for 3A; + 1 even : <^ t„, , -,. , ,, ^^^ U, , -,\ 

A3k+1) ^ >_(3fc4.3)/2 v^oo ~{3k+l) ^^s{i-l) 
' A3k+1) _ >._(3fe+3)/2 v^oo -(3A:+1) >_3(i_i) 

for 3A; + 1 odd : <; ^^,^^^ ^ ^.(3.+6)/2 ^p| gWi)^-3(.-i)! (^'^S) 

(C) ~ C^/^ V cf +\-i-3(.-i) ^ ^(3.+2) ^.^(2(3/2) + f(f^^^ cos(2C^/^), (5.36) 






for 3A; + 2 even 



' A3k+2) ^ >_(3fc+7)/2 y^oo -(3A:+2) >_3(i_i) 
^(3fc+l) _ >_(3A;+4)/2 y^oo g(3fe+2)^_3(j_i) 



{r(3k+2) ^ >-(3A:+4)/2 V^oo ^(3fe+2)^_3(j_i) 
f(3k+2) ^ ^-(3.+7)/2 ^pl gW2)^_3(.-l): (^-37) 

oo 

i^3.+3(C) ~ C^/^ E cS^'^^-^-^^^-^) + /f +^ sin(2C^/^) + fc'^' cos(2C^/^), (5.38) 

{p(3A:+3) _ >-(3fc+5)/2 V^oo ^{3fc+3) >_3(j_i) 
J S ~ ^ Z^j=l "-i ^ ' 

f (3A:+3) _ ^-(3fc+8)/2 y^oo g(3fc+3) >_3(j_i) 

{r(3A:+3) _ A-(3fc+8)/2 V^oo ^{3fc+3) >_3(j_i) 
^(3fc+3) ^ ^-(3fc+5)/2 y^'^ g{3fc+3)^_3(i-l)_ ^^'^^^ 

Moreover, suppose we have the following formula 



3 

AN 



-^0 i=0 



where e^ are all independent of A^ and 



(5.40) 



(5.41) 
Here the differential operator 

i (^) 

appears (3j — 1) times in the nested set of derivatives appearing in the integral above. 
According to the asymptotic expansion for Fsj{Q, it is easily seen that F^j^Q — c[ ^'(^^1'^ 
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is uniforinly bounded. This, together with (5.5), means that the error term in (5.40) 
satisfies |A-+i| <CN-^^^H). 



Now let us apply integration by parts again to derive next term in (5.40). According 



to definitions of -Ffc(C) in (5.32), it is easy to see that 

1 d 



^3,(0 -crc 



(3jVi/2 



C'dz 



Fsj+iiC), 



(5.43) 



Then integration by parts gives us 

""*+" d n d n d fe'(z) 
Tz [j'Tz VC 



A 



i+1 



.i)i+iA r 



dz V C 



"^zj+idz, 



(5.44) 



where now the differential operator (5.42) appears 3j times. Two more integration by 
parts gives us 



A, 



+1 



-1 



' AN 



z*+e 



d f 1 d 



dz \C'dz \C' 



d fe\z) 

d~z\Y 



Fsj+sdz, 



(5.45) 



where now the differential operator (5.42) appears (3j + 2) times in the above integral. 
According to Lemma [TI one can verify that the expansions in (5.34), (5.36) and (5.38) 



are also valid for k = j. Therefore, using (5.38) for k = j, we split (5.45) into the 
following two terms 



where 
e^.+, = (_i)J-+2^2,+i 

and 
A 



'"+' d n d n d (d'{z) 

Tz \X'~d~z \X''Tz \^ 



(5.46) 



^{3i+3)^l/2)^^ (5.47) 



■i+2 



-1) 



i+2. 



AN 



Z +£ 



d / 1 d 



dz \Cdz \C' 



d fe'{z) 

d~z\Y 



(i^3,+3 - cf^^''e/')dz. 

(5.48) 
Here the differential operator (5.42) appears (3j + 2) times in each of the above integrals. 



Recalling the definition of C, in (5.5 ) again, one can see that Cj+i is a constant independent 
of A^. Moreover, we have \Aj+2\ < CN'^^^^'al 



Therefore, we have shown that (5.40) also holds when we increase j by 1. As a 
consequence, (5.4) has an asymptotic expansion in powers of 1/A^^ and our lemma is 
proved. D 

Now we can prove our main results. 
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Proof of Theorem^ Let us recall (5.2). As we mentioned earlier, based on the results 
in [8] , the second term of its right-hand side satisfies an asymptotic expansion in powers 
of l/N"^, whose coefficients are analytic functions of t. According to Proposition p^ and 
Lemma |6| the first term of the right-hand side also satisfies an asymptotic expansion in 
powers of 1/N'^. The analytic property of coefficients in the expansion follows from the 
similar analysis as in [8]. Adding these two terms together, we prove our theorem. D 



Remark 3. As mentioned at the beginning of this paper, a combination of (1.10) and 
Theorem |2] gives us Theorem [TJ 
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